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Abstract 

We present an analytical derivation of the switching field distribution for a single domain particle 
from the Neel-Brown model in the presence of a linearly swept magnetic field and influenced by 
thermal fluctuations. We show that the switching field distribution corresponds to a probability 
density function and can be obtained by solving a master equation for the not-switching prob- 
ability together with the transition rate for the magnetization according to the Arrhenius-Neel 
Law. By calculating the first and second moments of the probability density function we succeed 
in modeling rate-dependent coercivity and the standard deviation of the coercive field. Comple- 
mentary to the analytical approach, we also present a Monte Carlo simulation for the switching of 
a macrospin, which allows us to account for the field dependence of the attempt frequency. The 
results show excellent agreement with results from a Langevin dynamics simulation and therefore 
point out the importance to include the relevant dependencies in the attempt frequency. However, 
we conclude that the Neel-Brown model fails to predict switching fields correctly for common field 
rates and material parameters used in magnetic recording from the loss of normalization of the 
probability density function. Investigating the transition regime between thermally assisted and 
dynamic switching will be of future interest regarding the development of new magnetic recording 
technologies. 
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I. INTRODUCTION 



The switching behavior of a single domain magnetic particle is strongly influenced by 
the presence of thermal fluctuations, which help the magnetization to overcome the energy 
barrier that is separating the two stable magnetization states. This leads to an effective 
reduction of the coercive fleld at flnite temperatures, which depends on the time scale of the 
experiment. The detailed knowledge of coercivity as a function of measurement time t and 
temperature is important for extracting the relevant magnetic properties (such as the zero 
temperature coercivity, i.e. the anisotropy H^, thermal stability ratio /3 = KV / {ksT)) out 
of experimental data. Sharrock pQ gave an analytical expression for a pulsed fleld experiment 



and others ([2], [3], [1], [5]) have also succeeded in deriving expressions for experiments in 
which an external fleld is swept at a constant rate. However, to our best knowledge, there 
is so far no explicit analytical derivation of the probability density function (PDF) describ- 
ing the switching fleld distribution (SFD) arising from the presence of thermal fluctuations 
exclusively. Knowing the SFD and its relevant parameters is of great importance for the op- 
timization of recording media as well as sensing technologies for magnetic flelds as described 
in [6]. 

According to the Neel-Brown model [7] , [8] the magnetization's spacial orientation fluctuates 
due to random magnetic flelds present at flnite temperature. This results in thermal insta- 
bility and forces the magnetization to switch between its two stable orientations separated 
by an energy barrier AE at a rate 



attempt frequency at which the magnetization tries to switch orientation and is usually 
treated as a constant. However, as already shown by Brown p], even for the very simple 
case of a single domain particle with a fleld H apphed parallel to its easy axis, the attempt 
frequency (in Hz) takes the form 




(1) 




(2) 



Eq. [2] is known as the Arrhenius-Neel law. The preexponential factor /q is called the 
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where a is the damping parameter from the Landau-Lifshitz-Gilbert equation, 7 = 
lei ii-^ = 2.21 ■ 10^ m/(A s) is the electron gyromagnetic ratio divided by the permittivity, 
Hx in amperes per meter is the anisotropy field, J5 in Tesla is the saturation magnetization, 
ks = 1.38 ■ 10^^'^ J/K is Boltzmann's constant and T in Kelvin denotes the temperature. 
In the first part, following the work of Kurkijarvi [Sj we will derive the switching field 
distribution from a master equation and the Arrhenius-Neel Law. In the second part, we 
will compare the analytical result to the output of a Monte-Carlo simulation written in 
MATLAB. The simulation allows us to introduce a function, which computes the attempt 
frequency according to eq. [3] in every time step. This way it is possible to compare the data 
from the Monte-Carlo simulation to a Langevin dynamics simulation based on the Landau- 
Lifshitz-Gilbert equation. 

Furthermore, we will discuss the rate dependence of the coercivity and its standard devi- 
ation and will compare the results to the models described in references [2] - [5|. In the 
conclusion we will address the open questions considering the validity of the Neel-Brown 
model, which are closely connected to the understanding of the transition regime between 
thermally activated and dynamic switching of a single domain particle. 

II. ANALYTICAL MODEL 

The energy landscape for the magnetization of a single domain particle with an external 
magnetic field H applied parallel to its easy axis has two stable minima separated by a 
barrier flU\ 



where K is the anisotropy constant and V the volume. There exist similar expressions for 
more complex reversal paths of the magnetization, which essentially differ in the exponent. 
The exact shape of the energy barrier can be calculated using for example the nugded elastic 
band method [TT] . 

When the field H(t) is ramped up, the time dependent probability Pnot that the particle has 
not switched until a certain moment to, is described by a master equation: 




(4) 



where f is the transition rate from the Arrhenius-Neel Law given in eq. [2j It has to be 
pointed out here that the Arrhenius-Neel Law only applies in the limit of AE > ksT as 
originally stated by Kramers [THj within transition state theory. 
From eq. [5] we get the following expression for P^, 
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as also derived in [5]|. To solve the integral on the right hand side of eq. [6] we use the 
substitution 

H{ty 



u 



K 



KV 



and 



Hit) = Rt 



where R is the field rate for a linearly swept field in T/s. The result of eq. Mis then 
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and the probability of switching is given by 
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R{uq{Hc)) is the cumulative distribution function (CDF) which describes the likelihood for 
the particle to have switched at a field H . The switching field distribution (SFD) is then 
given by the probability density function (PDF) which is the derivative of eq. [s] 
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III. RESULTS 



A. Switching field distribution 



In fig. |III A| the analytical SFD from eq. |9]is plotted together with the results of a Monte- 
Carlo simulation written in MATLAB. We assume a single macrospin switching between 
two magnetization states (up= +1, down = -1) in an external field H{t) swept at a constant 
rate. When the external field is varied as a function of time the switching probability is 
calculated by [TT] 



R = exp < — At/o exp 
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where At denotes the time step of the simulation and H{ti) is the value of the swept field 
at a certain time t,. Following the Metropolis algorithm, the values of P are then compared 
to a random number x G [0; 1]. The value of H gets accepted as a switching field He as 
soon as X > P. At a value H{t) > Hk the magnetization is automatically switched. We 
performed 100 sweep cycles and then computed the histogram of the obtained values for the 
switching field. The material parameters applied in the simulation are typical for modern 
perpendicular magnetic recording materials, /o was calculated using eq. [3] for zero field. 



We can see, that the simulation data agree very well with the analytical result from eq. [9] 
for the lower field rates. At high rates, however, a large discrepancy occurs between the 
analytical result and the simulation. If we compute the integral over the analytical SFD for 
8 T/ns we get a value < 1, which contradicts the general property of a PDF that 

""^ dP 

This implies that the CDF defined by eq. [8] must only give values between at — oo and 1 
at +00. For erf(+oo) = 1 we get 
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For the parameter system discussed here the right hand side of eq. H gives a value of 
7.2 T/ns, which explains why eq. M fails as a model for the simulation data at a rate of 



8 T/ns as shown in fig. Ill A 



In fig. Ill A we present results from a micromagnetic simulation solving the Landau-Lifshitz- 
Gilbert-equation with an additional random thermal field [I^ together with a modified 
Monte Carlo simulation, where the attempt frequency /o is now calculated in every field 
step according to eq. |3} To keep computation times resonable 80 mT/ns was the lowest rate 
that could be simulated within the micromagnetic approach. We see good agreement of the 
two simulation data sets (blue), however, fitting the analytical expression for the SFD from 



Chantrell et al. ^ R{Hc) = Hk/o exp[-/3(l - Hc/HK?]/[2f3{l - Hc/Hk)] 

Feng and Visscher [5j R{Hc) = {2\n2P)-^TT f^Hx {l - erf[V?(l - Hc/Hk)]] 

El Hilo et al. [3] Hc{R) = Hk{1 - \//3-i ln[/o/7x/(2i?/3)]} 

Peng and Richter 0] //c(iefr) = Hk{1- ln[/oteff/(2 In 2)]} teff = ieff(i2, /3) 



TABLE I. Analytical expressions for rate-dependant coercivity 

eq. [9] does not work satisfyingly. Fitting the data from the Monte Carlo simulation works 
better for lower rates, which is indicated by the decreasing standard deviation of the fit 
parameter /o,efr, which we call the effective attempt frequency. Compared to the model in 



fig. Ill A where /o was set constant to foi^H = OT) = 21.1 GHz we see a shift of the mean 
switching field to higher values, which is related to the lower effective attempt frequency we 
get from the fit for which fo^H = He) < /o,eff < fo{H = 0) applies. 



B. Rate dependence of the coercivity 

The coercivity of a magnetic recording medium is of special interest regarding the write 
process of a magnetic bit. Coercivity measurement methods usually employ very low field 
rates (e.g. VSM: <1 T/s) in contrast to the high field rates of up to 1 T/ns used in modern 
hard disk drives for the write process. Various models have been proposed to extrapolate 
between the laboratory time scale and the actual time scale used in magnetic recording. An 
overview of the expressions derived in references [2] - [H] is given in table |l] 

We calculate the switching field directly from the SFD. The mean value H of the variable 
H with the PDF dP/dH is given by 

/""^ dP 

Because dP/dH is not an analytical function we have to use numerical integration to calcu- 
late H. As it is shown in fig. [2] the switching fields derived this way are in good agreement 
with the other models, which have proven to adequately describe realistic recording media 
up to a certain field rate. We observe, however, that the models strongly deviate from each 
other at approx. 1 T/ns, which is a typical rate used in magnetic recording. As already 
stated above, the limit for failure of our model is determined by the loss of normalization 
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of the PDF (see eq. [TT). Peng and Richter |1] give R « /oi/x(21n2v^)-i = 2.9 T/ns 
as a condition for the vahdity of their model, which is a similar expression to the limit we 
derived in eq. 11 The expressions by El Hilo et al. P] as well as Peng and Richter [1] give 



non-real values as soon as the rates are out of the range of validity for the assumptions 
made to derive their expressions. 



C. Rate dependence of the standard deviation 

The standard deviation of the switching field can be calculated from the second moment 
of the PDF, i.e. the variance: 



vai{H) = a=J {H-Hf—dH (13) 
y J-oo dH 

Fig. [3] shows the standard deviations ac normalized to the switching fields Hq plotted 
with respect to the field rate. It is interesting to see that ac/Hc first decreases slightly, but 
then increases at higher field rates. Most importantly, ac/Hc has a minimum finite value 
at an intermediate field rate. We like to emphasize that the standard deviation in eq. [13] is 
derived for thermal fluctuations acting on a single domain particle, which would correspond 
for example to a single grain in a perpendicular recording medium. This is an entirely 
different effect as the standard deviation described by Hovorka et al. [14J , which arises from 
anisotropy and volume distributions of the magnetic medium. Hence, a direct comparison 
of the values for ac/Hc to values from the relevant literature is not reasonable. We are 
convinced, however, that the standard deviation derived here represents a fundamental limit 
regarding the accurate switching of nanoscale magnetic particles. 

IV. CONCLUSION 

We described an entirely analytical approach to derive the SFD of a single-domain particle 
caused by thermal fluctuations. Unlike the SFDs derived from an underlying distribution 
of grain sizes and switching volumes, the distribution presented here originates from the 
thermal activation of the particle's magnetization for overcoming an energy barrier at a 
rate described by the Arrhenius-Neel Law. The effect exists independently of any other 



contribution to the SFD and therefore represents a fundamental hmit to the accuracy of 
magnetization reversal at finite temperatures. An experimental evidence of the validity of 
the Neel-Brown model at low field rates was given by Wernsdorfer et al. [TB]. For typical 
field rates of > 0.1 T/ns used in magnetic recording micromagnetic simulation tools based on 
solving the Langevin equation describing dynamic switching under the influence of thermal 
fluctuations are used to compute the switching fields. However, computation time increases 
dramatically towards lower field rates, where thermal switching is dominant. It will be 
subject to future work to investigate this transition regime and understand the underlying 
physical processes. 

Furthermore, we are able to give an upper limit for the field rates where the assumptions 
of the Neel-Brown model are applicable. As with increasing field rates the switching fields 
approach the value of zero temperature coercivity, we show failure of the model due to the 
loss of normalization of the PDF describing the SFD. Using a Monte Carlo simulation for 
a single macrospin we reproduced the result of the analytical model within its range of 
validity and we were also able to include the field-dependence of the attempt frequency for 
the transition rate. By doing so, we observe significant changes to the switching field. We 
conclude, that the usual approach of taking the attempt frequency as a constant does not 
adequately describe rate-dependent coercivity, which is additionally supported by Langevin 
dynamics simulation data. 
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(a) Simulation data and analytical result for the switching field distribution, /o = const. 
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(b)Monte-Carlo simulation with /q computed according to eq. |3]in every field step and 
results from a Langevin dynamics simulation (a = 0.02). The solid lines are fits using 

eq. [9] with /o as the fit parameter. 

FIG. 1. Switching field distributions at different field rates for K = 0.3 MJ/m^, V = (8.8nm)^, 
Js = 0.5 T, fxoHK = 1.508 T 
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FIG. 2. Comparison of the models from table |l] and eq. 12 for the rate-dependent coercivity 
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FIG. 3. Rate dependence of the standard deviation of the switching field 
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